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NOTES ON TEE LAGRANGIAN MULTIPLIER METHOD 
IN ELASTIC-STABILITY ANALYSIS 
By Bernard. Budiansky, Pal C . Hu, and Robert W . Connor 


SUMMARY 


Nev applications of the Lagrangian multiplier method to stability 
analysis are described by means of elementary examples. The use of 
the method In analyzing the stability of (a) clamped plateB in shear 
cmH ( b) plate-stiffener combinations Is demonstrated. A detailed 
analysis for finding upper and lower limits to critical stresses of 
clamped rectangular plates Is presented in an appendix. 


INTRODUCTION 


The use of the Lagrangian multiplier method to calculate upper and 
lower limits to the critical compressive stress of a clamped plate was 
presented in reference 1. The procedures of reference 1 have been 
directly used to analyze the stability of clamped plates under com- 
pression in two directions (reference 2, footnote) and may be used, 
with little modification, to find vibration frequencies of clamped 
plates. The purpose of the present paper Is to describe additional 
applications of the Lagrangian multiplier method to the elastic- 
stability analysis of (a) clamped plates in shear and (b) plate- 
stiffener combinations. Elementary examples are used to bring out the 
essential features involved in applying the method to these types of 
problems. A detailed analysis for finding upper and lower limits to 
critical shear stresses of finite clamped plates Is given in an 
appendix. 


SYMBOLS 


a 

b 

X 

P 


length of plate 
width of plate 
half-wave length 
a/b or X/b 
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t 

E 

U 

D 


(^stiff 


A 


7 


El 

w 



plate thickness 
\ 

Young’s modulus of elasticity 
Poisson’s ratio 

JJ4-3 

plate stiffness in bending 

12(1 - ,* 2 ) 

effective flexural rigidity of stiffener attached 
to plate 

stiffener cross-sectional area 


a 

T 

k 


critical, compressive stress 
critical shear stress 

critical compressive stress coefficient in the 
formula a « 

critical shear stress coefficient in the 
formula t * k s 




x 

y 

w 


a.b.d, A 
n n' n 


*W b : 




n 

TKl 




plate coordinate parallel to length 
plate coordinate parallel to width 
deflection normal to plane of the plate 

Fourier coefficients 


Lagr an glan multipliers 


V 


internal bending energy 
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5 . 


J mn 


external work of applied. stress 
Integers 

Kronecker delta: 1 If m * nj 0 If 


A n . (l + l^n 2 0 2 ') 2 (l + 5^) 


B. 




n 


11 ^ " D n 2 

D n * UkgB^ 


En 


(l + n 2 S^) 2 - e 2 ^ (l + Son) 


2(7 - 3 2 Bk) 

Ann - 2^a 2 + n% 2 ) 2 (l + S^ + 5^ 


B 


v inn 


11111 V 2 - D im 2 




“ V 2 - D mn 2 


D nm * 


A'nm " (* 2 + ^23 2 ) 2 (1 + 5^) 


B* 


A* 


mn 


mn 


( A v) 2 - ( B, m) 2 


m n 



h 
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D, mn M ^kglB^mn 


Subscripts j 

pi plats 

stiff stiffener 

BUCKLING IN SHIAH OF AN INFINITELY LONG CLAKPKD PLATES 


Consider an infinitely long clasped plate loaded In shear as shovn 
in the accompanying sketch: 7 


The problem of finding the critical shear stress, for which tnere exists 
an exact solution by the differential equation approach (reference 3), 
trill be solved by the Lagranglan multiplier method. 

Boundary conditions .- It is evident that only the boundary condi- 
tions along the long edges affect the buokling stress. The boundary 
conditions of the problem are then: 

Zero deflection, long edges 


(//////;//////////////////, t/77777I7T//////f///////////////////, 



V/// ? ////////7/////////;//77T ////////////77/////7/////////77, 




Zero slope, long edges 


( 1 ) 
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Fourier expansion .— The deflection surface of the buckled plate Is 
known to be harmonic along longitudinal sections; the deflection surface 
may therefore be represented in the form 


v = fi(y) 8 in y + f 2 (y) 00s y (3) 


where \ is the half-wave .length. The necessity remains of choosing 
suitable Fourier expressions for f -j_( y) and fg(y) for use in the 

Lagrangian multiplier method. It is desired that the series chosen 
satisfy the conditions that: 

(a) In the region 
symmetrical (see preceding sketch) about -tie origin 

(b) The potential energy expression for the buckled plate, calcu- 
lated on the basis of the specified expression for w, consist of 
integrals of products of functions which form orthogonal sets. 

An expression satisfying these requirements is 

w = sin y ^ Oq sin 

n*l n=»0 


2nny 

b 


+ cos 


CP 

fZ 


djj cos 


2n«y 

b 


W 


V 2 * 2 .)* 


the deflection w be 


The fulfillment of condition (b) is verified in the following section, 
in which the energy expressions are calculated. The desirability of the 
condition will become evident when the final stability criterion is 
derived . 

Energy expressions .— Substituting the value for w from equation ( 4 ) 
into the formulas for bending energy and work 





- 2(1 - u) 


5^w 8^w 
3x 2 8y^ 



dx dy 


( 5 ) 


T = 


— Tt 


dx'dy 


( 6 ) 
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Constraining relationships .— Neither of the clamped-edge boundary 
conditions (l) and (2) are satisfied term by term by the expansion (4)j 
constraining conditions must therefore be imposed on the Fourier coef- 
ficients a n and d^. These constraining relationships are, for zero 

deflection, 

CP 

(-l) n dj, - 0 (8) 

n-0 1 

and for zero slope. 




n*l 


(9) 


Stability criterion .— The energy method requires that V — T be 
minimized with respect to the a's and d*s. Since the a*B and d's 
are, however, bound by equations ( 8) and (9)> the minimization is 
performed by the Lagrangian multiplier method by minimizing 


UJ 

■ ^ IT (-l) n dn - t n(-l) n a T1 (10) 

( ^ 1 ^ ^ — 

Wv 


n*0 


n»l 


with respect to the a* s and d's. The Lagrangian multipliers are p 
and £. The complete set of equations for minimizing V — T with 
the a's and d's bound by the conditions given in equations (8) 
and ( 9 ) then becomes 


- 0 (n - 1,2,3,...) (11) 

8an 



- 0 (n * 0,1,2, . . .) 

CD 

(-D n ^ - 0 

n=0 


n(-l) n a n = 0 

n=l 


( 12 ) 


( 8 ) 


(9) 
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Substituting aquation ( 7 ) into aquation (10) and aquation (10) into 
aquations (ll) and (12) gives 


2(1 + W a n 2 ) 2 a n - C (n - 1,2,3,...) ( 13 ) 

- (SV 3 ”)*!! + + So^dn . (-l) n ,, (11) 

(n - 0,1,2,...) 


Solving aquations (13) and. (14) simultaneously for a n wtvi gives 


*4i »§[B n n(-i) n C + c n (-i) n u] 
*4i “ 2 [pn n (“*^-) n C + 


(15) 

( 16 ) 


where 


Bn 


An 2 ~ D n 2 


11 "Bn 2 


in which 


An - C 1 + 4n2 P 2 )(l + 0 


D n » fckgp3n 
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Substituting equations (15) and (l6) into equations (9) and (8), 
respectively, yields 


p S~ ^TL + nC n - 0 

n«0 n«l 




^n 



n*l 


0 



(17) 


The condition that a nonvanishing solution exists for |i and £ gives 
as the final stability criterion 



For a given value of P, the value of k a that satisfies 
equation (18) can be found by trial substitution and interpolation. 

The correct value of P is that which gives the lowest value of k 0 . 

For P = 0.80, equation ( 18 ) yields kg = 8.989 which agrees with the 
solution obtained by Southwell by the differential-equation approach 
(reference 3). 

Discussion of method; further applications .— In the usual appli- 
cation of the Rayleigh-Sitz method, an infinite set of equations 
involving infinitely many deflection coefficients is obtained when the 
energy expression V — T is mini mi zed ; the exact stability criterion 
is then an infinite determinant obtained from these equations . The 
simplicity of the solution Just obtained, however, is due to the fact 
t ha t it was practicable to transfer consideration of infinitely many 
Fourier coefficients to consideration of only two Lagrangian 
multipliers. An essential feature of the solution that permits this 
simplification is the fact that the substitution of the expansion 
chosen for w (equation ( 4 )) into the expressions for V and T 
(equations (5) and (6)) leads to integrals of products of functions 
which form orthogonal sets. It then becomes a simple matter to solve 
the minimizing equations ( 13 ) and (lU-) for the Fourier coefficients In 
terms of the Lagrangian multiplier (equations ( 15 ) and (16)), substitute 
back into the constraining relationships (8) and (9), and derive the 
stability criterion from the condition that there be a nonvanishing 
solution for the Lagrangian multipliers. 
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Stability problems involving finite rectangular plates may require 
double rather than single Fourier expansions for use in the Lagrangian 
multiplier method . However, a proper choice of series, with due 
attention to orthogonal! ty considerations, would still make It possible 
to solve for infinitely many Fourier coefficients in terms of 
Lagrangian multipliers. Appendix A presents a solution by the 
Lagrangian multiplier method for the problem of the buckling in shear 
of a rectangular plate clamped' on all sides.' Although it was possible 
to obtain an exact solution for the infinitely long plate, the solution 
for finite plates is approximate; however, as in the compress! ve 
buckling problem discussed in reference 1, it is possible to obtain 
approximate solutions in two different ways, which permit the compu- 
tation of an upper limit as well as a lower limit to the true critical 
stress. The true critical stress can thus be bracketed to within any 
desired degree of accuracy by taking sufficiently close upper— and 
lower— limit approximations. 

The Ijagrangian multiplier method may, with the use of appropriate 
deflection functions, find applications to other problems. A general 
discussion of Fourier series and their use in stability analysis is , 
contained in appendix B. 


BUCKLING OF PLATE-STIFFENER COMBINATIONS 


The application, of the Rayleigh-Ritz energy method to buckling 
problems involving plates with stiffeners usual. ly results in energy 
expressions that are conpllcated functions of the deflection coeffi- 
cients. That is, even if the terms of the assumed deflection function 
have the orthogonality properties previously discussed, energy terms 
due to stiffener deformations will usually involve quadratic cross 
products of all combinations of the deflection coefficients. 
Occasionally, for some special problems (see reference If) relatively 
simple stability criteria can still be derived by algebraic manipula- 
tions; however, in general, it is to be expected that an exact sta- 
bility criterion for a stiffened plate, derived by the Rayleigh-JRl tz 
method, will consist of an infinite determinant that is obtained from 
explicit consideration of infinitely many deflection coefficients. 

The Lagrangian multiplier method can be used to simplify the 
analysis considerably. As in the unstiffened plate buckling problems 
previously discussed, an appropriate application of the Lagrangian 
multiplier method makes it possible to Bolve for Fourier coefficients 
in terms of Lagrangian multipliers, so that explicit consideration of 
a finite number of Lagrangian multipliers takes into account infinitely 
many Fourier coefficients. The elements of the method of application 
of the Lagrangian multiplier method to stiffened plates will be- pre- 
sented by giving the analysis of the stability under longitudinal 
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“V 


compression of an infinitely long damped plate haring a longitudinal 
stiffener alon g its center line as illustrated in the following sketch: 


/ 

///////////, 

\ Stiffener ~7 

Zzj///////////////////y////////////////// / 

A 

* 



4 

/ V 

•< X 5- 

}))')))))))))) 

7-m-rrf 


In this example, the stiffener is assumed to hare no torsio nal rigidity. 
(See reference 5 For solution of this problem by the differential- 
equation approach.) 

Boundary and continuity conditions .— The boundary conditions along 
the clamped edges are zero deflection, 

v pl( x *°) “ v pl( x > b ) -° (19) 


and zero slope, 

^( x ,0) » ^U,b) * 0 (20) 


The condition that there be continuity between the plate and 
stiffener is given by 

w stiff “ v pl * 0 (21) 


Fourier expansions The buckled deflection surface is known to be 
sinusoidal in the long direction. However, the deflections in the short 
direction may be either symmetrical or anti sysmetri cal, depending on 
which mode corresponds to a lover buckling stress. For the present, 
the symmetrical mode will be considered. Then, let 

Vpi = sin 2 S hfc cos ajj* 

n=0,2,4, . . . 


( 22 ) 
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and 


v stiff 


A a in ^ 


(23) 


where X La the half-wave length of the 'buck: led surface. 

Energy express Iona In addition to the plate— bending energy 
(equation 5 ) , the stiffener-bending energy 


(El) 


stiff 


/ d a v, 


r stiff 


stiff 
dx 2 J 


dx 


(24) 


the work of the plate stresses 


T 


Pi 



dx dy 


(25) 


and the work of the stiffener stresses 


1 „ OA 

stiff ~2 



dw 


stiff 
dx J 


dx 


( 26 ) 


must he taken into account. Substituting equation (22) into 
equations (5) and ( 25 ) and equation (23) into equations (24) and ( 26 ) 
gives, in the region (0,\),(0,b), 


v pl + v stiff “ T pl “ T stiff 


itS)b 

8X 3 



n-0,2,4,. 


(l + n 2 p 2 ) 2 - e 2 * (1 + So^bn 2 


+ 2(7 - e 2 5k) A 2 


(27) 
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where 



Constraining relational 1 ^ pa .— Will la the 'boundary conditions of 
zero slope (equation (20)) are satisfied by each term of equation (22), 
the boundary conditions of zero deflection (equation (19)) are satisfied 
only by making 

CO 

ZL bn - ° (28) 

n*0,2,lv, • • • 


The continuity condition (equation ( 21 )) will be satisfied by means of 
the constraining relationship 


A - 1 ^ cos 2 ^ « 0 

n* 0 , 2 , 4 , ... 

or 

A - ^ (_i) n / 2 b n « 0 

n*0,2,k, . . . 


( 29 ) 


Stability criterion .— The energy expression (27) must be minimized, 
with A and the b’s bound by the constraining relationships (28) 
and ( 29 ). Let 
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where § and. T) are Lagrangian Multipliers . Then, 

g- - 0 - 2^(l + n 2 0 2 ) 2 - P 2 * (l + S 0n )b n - 5 + n(-D 
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Substituting aquations ( 30) and. ( 31 ) into equations ( 28 ) and (29) gives 


E — T1 
n ' 


n- 0 , 2 ,^, . . . n=0,2.,k, . . 


Y_ (- 1 ) 11 E n * 0 


y (32) 


- 5 2 ^ (-l) n 'Zn + n (h + ®n) = 0 

n*O t 2,k, . . . V n= 0 , 2 ',^, . . . 


Setting the determinant of the coefficients of I and r\ in 
equations ( 32 ) equal to zero gives as the stability criterion for 
symmetrical buckling 


!>_ ** 

H + ** 

_ 

l 

M* 

r, 

11=0 J 2 f ^ y m m • 

n=0 J 2,l)- J . . . 


n=0,2,k, . . . 


(33) 


For a given value of p, this criterion can be used to find kj the 
correct value of P is that which gives the lowest value of k. 

The stability criterion for anti symmetrical buckling is obtained 
simply by letting 

m 

v * 8in X 21 COB X 1 

n=l ,3 » 5 » • • • 


and using the Lagrangian multiplier method to introduce the zero 
deflection constraining relationship 



n- 1 , 3>5j • •• 
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(Since the stiffener has no torsional- rigidity and lies along a node, 
it does net affect the "buckling stress.) The stability criterion for 
antisymmetrical "buckling becomes 



n=l, 3 ) * • ■ 


^n - 0 


( 3*0 


If equation (3*0 gives a minimum value of k that is lower than that 
obtained from equation ( 33 ), the stiffener will remain straight and 
the plate will buckle- antisymmetrically. 

Further applications . - It is evident - that the method given may "be 
used to take into account any number of stiffeners with arbitrary 
spacing. Furthermore, the affect of stiffener torsional restraint can 
also be included by adding to - the energy expression the twisting energy 
of - the stiffener and introducing conditions of rotational continuity 
between plate and stiffener. The method can be extended to analyze 
stiffened plate problems in which double Fourier series must be used 
(for example , a finite rectangular stiffened plate, clamped along all- 
edges )j in such problems, approximate upper- and lover -limit solutions, 
rather than exact-solutions, may be expected. Stability problems 
Involving p late-s t if fener combinations in shear, or in combined 
compression and shear, might also be conveniently solved by the present 
method. 


CONCLUDING REMARKS 


Elementary examples have been given to demonstrate the application 
of the Lagranglan multiplier method to the elastic -stability analysis 
of (a) flat rectangular clamped plates in shear and (b) plate -stiffener 
combinations. Exact solutions were obtained for the examples con- 
sidered} for other problems, such as the shear buckling of a finite 
clamped plate (appendix A), approximate solutions may be obtained in 
two different ways providing upper and lower limits to the true value 
of the buckling stresses . 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., September 12, 19*+7 
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APPENDIX A 


BUCKLING IN SHEAR OF A CLAMPED RECTANGUIAR PLATE 


A thin rectangular plate clamped along al 1 four' sides and loaded 
In shear along the edges is illustrated in the following sketch: 


y 



The problem is to determine the critical value of the shear stress 
under which a plate of given aspect ratio begins to buckle. 

Approximate analyses of this problem have been given by Smith 
(reference 6) and Iguchi (reference 7 ). Smith uses the Rayleighr-RItz 
method and hence obtains an upper— limit solution, whereas Iguchi uses 
a method that leaves the type of solution unspecified . Both upper 
and lower limits to the true buckling stress may be obtained by the 
present Lagrangian multiplier solution. 

The buckling configuration may be either symmetrical or anti- 
symmetrical about the plate midpoint, depending on which buckling mode 
corresponds to a lower critical stress; the two buckling patterns will 
be considered separately. 


Symmetrical Buckling 

Boundar y co ndition s.— The boundary conditions of the problem (see 
preceding sketch) are: 

Zero deflection, all edges: 


w(x,|) -v(x,-|) 


0 


(Ala) 
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w(|,yj 

( = w ( _ | >y) ■ 0 

(Alb) 

Zero slope, all edges: 





(A 2 a) 

s(l ’ J ) 

■ s(*l ,y ) ■ 0 

( A 2 b) 


Fourier expansions In appendix B, conditions to "be considered in 
choosing Fourier series for use in stability analysis are discussed . 

On the basis of this discussion the following expansion was chosen to 
represent the symmetrical buckled surface : 

CP 00 

\ ^sT" , 2nutr , 2nity 

» - z_ 2 — 8ln — sr Bin nr - 

m=»l n=»l 


2i- 

m=»0 n=»0 


cos 


2mnx 


COB 


2n?ty 


(A 3 ) 


Energy expressions Substituting the expansion for w (equation (A3)) 
into equations (5) and ( 6 ) gives 



m=l n=l 
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Than 



V“ 2 (^ + n%2) 2 (l + 5^ + Sq*) 

Not© that V — T is independent of cLqq since Aq 0 - 0. 

Constraining relationships In order to satisfy the "boundary con- 
ditions of zero deflection (equations (Ala) and (Alb)) it is necessary 
to impose the constraining relationships 



m«0 


(-l) m d 


mj 


0 


(J - 0 , 1 , 2 ,...) 


(A5a) 


(-l) n d-m - 0 (i - 0,1,2,...) (A5b) 

n»0 


S imilar ly, in. order to satisfy, the zero slope conditions (equations (A2a) 
and (A2h)), it must he true that 


/^~ “( - 1 ) 111 V ) - 0 

m»l 

0> 

y n( - l)n a in = 0 
n=l 


( J =» 1> 2, 3,» • • • ) 


( i = 1*2, 3j • • • ) 


(A6a) 


(A6h) 



20 


NACA TN No. 15 5^ 


Equations (A5a) and ( A5b) contain the term clgg, which Is missing 
from the expression for V — T. The term dgg may he eliminated from 

the constraint relationships by subtracting the first of equations (A5b), 
the equation for i = 0, from the first of equations (A5a), the equation 
for J * 0. The final set of necessary constraining relationships then 
becomes 


S' <-D” <w> - Si ( - i)n ^ - ° 

m=*l n«l 


/ (-1)" - 0 

m*0 


CD 



n«0 


(-1) Q *in 


S= 0 


(J « 1,2,3,...) 
(i = 1,2,3,...) 


y (A7a) 


^ m (-l) m j “ 0 

m=l 

00 

^n(-l) n a ln “ 0 

n«l 


(J =■ 1,2,3,-..) 


( 1 = 1 , 2 , 3 , » • • ) 


(A7b) 


J 


Lower— limit solution .— As described in reference 1, a lower lim! t 
to the buckling stress can be found by minimizing the energy 
expression V — T (equation (Ah)) with respect to all the coefficients 
and < d inii but satisfying only some of the constraining relationships (A 7a) 
and ( A7b) , say, as far as i = p and J = q. 


If V - T 


is di Tided through by the constant term 


a 3 ' 


the 


function to be minimized becomes 
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v 1 ' S m0 ~ 5 0n, 


- - Z ( - 1)m - Z ( - i)n< 


-a. co -er 00 

Z x j Z ( - i)m Vi - Z “i Z £ - i)n a i 

J=1 m=0 1=1 n=0 


SL 00 ir co 

' Y- ^ Y- m( “ 1)m V) - I] nf " 1)I1 


j=l m=l 


i=l n=l 


The quantities a, \y u and are Lagrangian 

multipliers. The minimizing equations "become 


(m,n = 1 , 2 , 3, • • • ) 




(m,n = 0,1,2, ... ) 


Equations (A7a) and (A 7 b) taken up to 1 = p and J = q 


By evaluation. 




“ - k s p3mnd rnn - ^n m ^ -1 ^ m - = 0 


— Anrn dmn ~ k^B^mnarrm - Xn( - l) m - Um(-l) 11 = 0 (m,n 0) (All) 


= WmO - *(-V m - Vm = 0 


(A 12 ) 
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55 “^ = AOn^n + <*(-l) n - ^ = 0 (A13) 


In equations (AlO) to (A 13 ), £ n and X n do not- appear for 
values of n > q, nor do r) m and appear for values of m > p. 

V/hen both m > p and n > q, £ n , ri m , and vanish from 

equations (AlO) and (All). Then one of two conditions is possible, 
either 

A^ 2 - h s 2 3 6 m 2 n 2 = 0 


or 


a, ,,, ., ~ d = 0 

^TITI “ Wl 


The first alternative, however, for given values of m >'p 
and n > q will ordinarily lead to a very high value of the buckling 
stress coefficient kg. For the lowest buckling load, therefore, 

when m > p and n > q. 


a,™ = d^ = 0 
Tnn tnm 


gives 


For the remaining a 's and d's, solving equations (AlO) and (All) 


a mn — ®mn 


m(-lf ?n + n ( _1 ) n T lm 


+ C 


mn 


(-if ^ + (-l) n Ma 


(Al4a) 


<3inn = Enm 


(-if ^ + (-if 


+ C nm 


m(-l) m £ n + n(— l) n x\ 


■m 


(Allb) 
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where 


Arm = 


Amn 


2 2 
Amn — Ann 


Cmn - 


Arm 


. 2 _ 2 

Arm Prim 


Arm = i# 3 ™! 



ITroiii equations (A12) and (AI 3 ), 


= ^|a(-l) m + ^ 
4on ' ^ 


(A15) 

(Al6)‘ 


It is to he emphasized that in equations (Alta) to (Al6), for 
values of m 5 p 


v ■ ‘V ■ 0 


and for values of n > q 


=n 


= ^n - 


Substituting the values of a 's and d 's given by equations (Alta) 
to (Al6) back into the constraining relationships (A7a) and (A7b) taken 
up to j = q and i = p gives 



a ( S' ho + ]t_ 3 0n) + X- ( ~ 1)mB ^“21 ( " 1)riB 0^n =0 


X j2Z\j a (“ 1 ) J B 0J + 1)1B+J B mj u m + ^J^L^mJ “0 (j = 1,2, 3,.. c 


q <p q 

Z ( - 1)1+ " B ln\i * a( - 1 > 1 B 10 + *1 Z B in + Z 1( 


(— in c. ltrn ''1 / "In 

n=0 n=l n=l 


■O 

-z 


- X) C ir£n + \/ *£.=0 (i* 1,2,3,. .P)U(A17) 


«o P re P 

£ jZ A «J + Z^’^ B mj\ + \)Z *X) + Z V. * ° (J ’ 


* w q w 

22in(-l) i+n B ln C n+ 1^21 n% ta + ]r n (-l) 1+n C ln X n + " 0 {l = 


• 

In order for this set of 2(p + q) + 1 equations to be compatible, the determinant of the coefficients 
of the lagrangian multipliers must vanish. This requirement leads to a detenninantal equation from vhich 
the critical value of the buckling coefficient may be found by trial. An example of an eleventh-order 
determinant, with p - 3 and q * 2, is shown in table 1. 
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Upper— limit solution .— The theory of the upper-limit solution in 
the Lagrangian Multiplier method (reference 1 ) requires that some a’s 
and d*s arbitrarily be set equal to zero, that expression (A 4 ) be 
minimized with respect to the remaining a* s and d's, and that all 
the constraining relationships (A7a) and ( A7b ) be satisfied. 

A^ a result of the necessity for satisfying all the constraint 
relationships, a redundancy exists among equations (A7a) (see 
reference 1) ; this redundancy can be removed by discarding the first 
of equations (AJa) . That another redundancy exists In equations (A7b) 
may be shown as follows: If 


*BlJ * 0 


m=l 


Is multiplied by j(— 1)^ and summed over j, the result 1 b 


w 

jV 


.m 


V 


j(-l V 2L H ~ X) a mj = ZL Z_ Jn( " 1) V) = 0 (A18) 

J=1 m=l J =1 m=l 


. J+m 


and if 


n(— 1) a in = 0 


n=l 


Is multiplied by i(— l) 1 and summed over i, the result Is 


00 CO 


ZL 1(_1)1 zL n(_1)I1 84n = X I 1 *- 1 ') 1 ™ a in = 0 (A19) 


1=1 


n=l 


1=1 n=l 


Equations (Al8) and (A 19 ) are identical; hence a redundancy exists, 
which may be removed by discarding one of equations (A 7 b), for example, 
the equation for 1=1. 

With the elimination of the redundant conditions, the necessary 
constraint relationships become 
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y~ (-D m - ° 

m=0 

OP 

/" (-D n <kn - 0 
n=0 


(J - 1,2,3,...) 


r ( A20a) 


(i « 1,2,3,...) 


00 

* 0 

m=l 



n=»l 


(J - 1 , 2 , 3 ,...) 
(i = 2,3,^,-..) 


Y (A20b) 


The accuracy of the upper-limit result, as veil as the ease of 
solution, depends in part on which Fourier coefficients are retained 
in the analysis . Several possible groupings of included terms were 
tried out in solutions for the special case of a square plate. The 
trials indicated that the optimum arrangement for practical applica- 
tions was a finite rectangular array of cosine coefficients (^mn) 

together irith certain infinite rows and columns of sine coeffi- 
cients (djnn), as illustrated in table 2. 

These limits on the existence of the coefficients can be expressed 
as follows: 


dm - 0 (when either m > p or n > q) 


Vi 


(when both m > p and n > q) 


When these limits are imposed, the constraint relationships (A20a) 
and ( A20b ) take the form 


zL ( ~ 1)m Sj 0 

m=*0 


(j * 1 , 2 , 3 , • «q) 


(_D n d ln - 0 (i = 1 , 2 , 3 ,.. p) 

n-0 


r 


(A21) 
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<33 



(j * l#2#3i • # l) 

m-1 

►* 

P 



( 0 ■ • • • ) 



^ n(-l) n a ln - 0 

(i » 2,3,h,..p) 

n=l 

r 

<1 


^ »(-l) n a in “ 0 

(i » p+l,p+2,...) 

n=l 



The function to he minimized. 

is 



OD CO 

-XX 

|"^rnn[®m 2 ( 1 

®n O ” 1 

- 


m-0 n-0 

L 



-J 

-±'*: 

^ Vj - 

P 

Z-i 

X t_i)n ^ 


J=1 m-0 

i-1 

n®0 



00 

-X 


Vl 

J=1 m=l 

3=*q+i 

m=l 


P CO 

^i y~ n(-D n ®i n 

i=l n=l 

03 

-X 

i =p+i 

q. 

Hi n ^“ 1 ) n 

n=l 

*in 


(A22a) 


( A22b) 


(A23) 


Sotting 


as 


ran 


5d 


mn 


0 then gives 
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For 4 jnn * °> 

«S=a " B «[“(- 1 ) B 5n + a (~V n %] + Gan^- 1 ) 111 *n + <-^> n Mm] 

For - 0, 

“no “ £(■(-*>“ «n + n(-l) n tJ 


Substituting these values bach into constraint relationships (A21) 
through (A22b) gives 
» 

P P P 

X J y + + 2Z 

m =0 m =-0 m >=0 


P_ 

+ ]> (-D nM ' J JC oJ Ti m = 0 (j = 1 , 2 , 3,.. l) 
m=0 


q. 

Z<- 

n=0 


•1) 


1+n 


B. X. 
in n 


+ 


H. 


1) 


i+n 


iC in*n 


n=0 


n=0 


+ 1i 


q 

X 

n=0 



a 0 


(i = 1,2,3,. -p) 


P " 2 P 


Z ” aB »j Z ( - 4) “ J •*-}'- 

m=p+l “« m=l 


m=l 


Z^% t x 

m=p+l m=l m=i 


+ x 3 ^ “CmJ + 2I ( - 1)DH ' J 


( A2Ua) 


(A2kb) 


(A25a) 


fj - 1.2,^...^) 
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» \ m 2 V" mj 


( J - q+l,q+2,. ..) 


£ <-« 1+n z 

11=1 n*q+l 



+ n. 


I £*X ( ~ : 


■1) 


l+n 


nC X + 
in n 


n=q+l 


n=l 


4 . 

iZ 

n=l 


nC 


in 


(i = 2,3,4, ..p) 



(i * p+l,p+2, .. .) 


These equations involve all the multipliers £ j and 

can he reduced, however, to a set of equations including only 
multipliers up to i = p, j = q as follows: 

From equation (A25b), for J = q+1, q+2, . . . 



From equation (A26b), for i = p+1, P+2,... 


^i 


q. 

-Z 

n=l 


(-D 1+n 

Ain 



n=l 



Ain 


-29 


(A25b) 


(A26a) 


(A26h) 


They 

the 



sm- 


Substituting th©a© ©xpreasions for and m Into equations (A26a) and (A25a) respectively 
.gives, along with aquations (A24a) and (A24b), the final stability aquations: 


1l Z n2a i» * Z^" V. + Z nl( - 1)1+n B la ? 


n + i*j } a - in + \ 


« 

7 n - 

A Air. 


n.nq+1 


P « 


min^(-l) nH ' 1 


m=2 n=q+l y~ 

A_ A_ 


= 0 (i = 2,3,4,. .p) 


! J Z A «J V» * l jZ A) + Z “ ( - :L) " J Vrn + l l Z 


m=p+l 


-ii 

n=l m=p+l 


• m 2 nj(-l) n ‘ f ^ 


Z_ A„ 


0 ( J = 1,2,3, . .<l) 


P P P 

X J Z_ Bm J + X_ (_1)mfJ B “j^ + + ZZ ^- 1 ^ J Cjnj’lm a 0 (J = 1,2,3,.. d) 


* i+n q q 

B in + 2_ (_1) +nB ln X n + ZI l( ^ 1)i+nC ln 5 n + ^if 1 ~ 8 li) ZZ “° {i aJL ' 2 *3,--P) 
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The stability criterion, is a determinant based on equations (A27), 
as in the lower limit solution. An example of a ninth-order deter- 
minant, with p * 3 and q. =« 2, is shown in table 3. 


Anti symmetrical Buckling 


For the analysis of buckling anti symmetrical about the plate 
midpoint, the origin of coordinates was taken as shown in the following 
sketch: 


7 


i 

k 



x 


Boundary conditions The boundary conditions are now expressed as: 
Zero deflection, all edges 


K 1 *!) - v ( I,_ s) * 0 

w(o, y) - W (a , y) = 0 

Zero slope, all edges: 

dw/C b\ dw/ bN 

SFVf'zJ * 3?(. x ' ~z) = 0 

ll ( °' t) ‘ ll (a * 71 “ 0 
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Fourier expansion . - In accordance with the Ideas s ummar l zed in 
appendix A, the chosen expansion is 


V 


00 00 

/* "^mn 
m=l, 3 a=0 



EBtX 

COS COS 

a 


2r aty 
h 


00 00 

+ XZ ZT Cmn 

m=l,3 n=l 


sin 


nutx 

a 


sin 


2mty 

b 


(A28) 


Energy expression .- When the expansion (A28) is substituted in the 
general formulas (5) and (6), the result for the total energy expression 
is 


Y 


T 


_ Pit;' 4 ! 
4a3 


00 00 

< /* ^ 5^ mnjZna 2 + c nm. 2 f 1 " 5 0n) 

- 

^m=l,3 n=0 L x 

J 


(A29) 


where ' 

- (» 2 + + 8 0n ) 


Constraint relationships . - The necessity o*f satisfying the 
boundary conditions ixrposes the following constraining relationships 
on the coefficients: 


00 

Y~ f- 1 ) 11 *in = 0 Ci = 1,3,5, •••) 
n=0 

CO 

Y~ *mj =0 (J = 0,1,2,...) 

5^173 


> (A30a) 
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y~ n(-l) n Cln . = 0 (i = 1,3,5,- ••) 
n=l 



m=l,3 


“mj 


= 0 


(j = 1,2,3,...) 


>> (A30b) 


Note that there la no term h 0 Q to he removed from these egressions. 
Lower— limit solution.— The function to he minimized becomes 


G- ' = 


w w 

n 

m-1,3 n=0 


2 *’ 

2^ Bin 


w 


+ c. 


mn 



P 

- tkgB^nh^c^ - J 

1=1,3 


op 

i YL ( - i)1 


j 1 xl 


n=0 


q. 

n 

J=0 


c 

I 


"Vj 


m=l,3 


) ti \ ^ n(-l) 31 
i=l,3 


'in 


n=l 


- 1 


me. 


mj 


J=1 m=l,3 



The stability equations analogous to equations ( AI 7 } for the symmetrical case talcs the form 


00 

'*X 


B» . + 


m=l,3 ®-l,3 


CO P 

»V’« 4 ^’jC 1 - 8 oj) X ,c Vj 4 X V* - 0 


m=l,3 m=*l,3 


(J = 0,1,2,.. q) 


q 00 q ** 

X<- 1)nB Vn * “’l X B ’ln + X l( - 1,n °'ln«* n 4 "’iX^'ln ’ 0 (1 ’ 1 - 3 - 5 -’ 


~ ±f W J- 

X m 2 B, »j + X b v » + x 'j X *'-) 4 X ° v» ■ ° 

m=l,3 m=d,3 m-d,3 ra«l,3 

( J * 1*2,3/ • «q) 

i 80 q. 

X (-1)“ inB' ln t' n 4 i,*i 2_n%« ln . X"*- 1 ' 11 + M’iX - 0 

n=l n=l n=l n«J. 

(i * 1,3,5,. .p) 


From these equations a determinant can be formed to give the critical value of the buckling 
coefflci ent. 


i 


4 &L VOVtt 
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Upper— limit solution .— Since all of constraint relations (A30a) 
and. (A30b) are to be satisfied according to the theory of an upper- 
limit solution, It can "be seen that one of equations (A30a) and one of 
equations (A30t>) are redundant . These redundancies can be proved in 
the same manner as for the symmetrical case and are eliminated by dis- 
carding the equation for i = 1 of equations (A30a) and also the 
equation for «J = 1 of equations (A30b). As was done in the upper- 
limit solution for the symmetrical case, a rectangular array of cosine 
coefficients and infinite strips of sine coefficients will be retained; 
thus, let 


and 


bjjjji = 0 (when either m > p or n > q) 


c mn “ 0 (when both m > p and, n > q) 


The final set of constraining relationships then becomes 



n=0 


(-D n v - 0 


z 

m=l,3 


V] “ ° 


n(-l) n * in - 0 


n*l 

4 




1 ) r, la ,o 


n*=l 


/ 

m=l,3 

m=l 


= 0 


= 0 


(i - 3,5,7,.. P) 
(j m 0,1,2, . .q) 


(i * 1,3,5,.. P) 


(i * p+2,p+h, ...) 

( J = 2,3>, • •<! ) 

(J = q+l,q+2, ...) 



Performing operations similar to those for the symmetrical case glyes as the final stability 
equations : 

q q q q 

Vi 2 >„ *zl B' ln t» n *y n (-i) n * li’ifi - C lt )/ °c'm 


m ¥ m 

n=q+l m=l,3 n=q+l \ m z 


n^mi . 

^ n ■ 


A' ia A 


1 — V 


> jf 

-C— A' 


m=l,3 


P P P P 

X X * l 'j Z " c '^ + 2_ 

m=l ,3 111=1,3 n=li 3 m= 3,5 


C* ,n' 


- P q oo » ^ JbU n _ m 

, V « V V A m/ V 5 n 

^ z. z_ "jl : =0 u “ 2,3A - ,q) 

m=p+2 n=2 m=p+2,p+4 \ ^ 

Z_aT~ 


r-i- a q q 

3 *in X ’n * 2 l I( ' 1)n C, ^'“ + ’’l H "°’ta ’ 0 !l ■ 3 ’ 5 ' 7 *-'' 


P P P P 

X ’j Y B Vi* Y (VuV.-tt'j^l-Ooj-Sij) X "C'»J* 21 J(-D J C'mjV = 0 

m=l, 3 n=3,5 H“l, 3 m=l,3 


( J - 0,1,2, ..q) 
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Numerical Example 

Upper and lower limits to the symmetrical buckling stress of a 
square plate have been computed; certain simplifi cations are possible 
in the operations for that case . From the symmetry about both diagonals 
present in the buckle pattern of a square plate, it can be shown that 


a mn = a nm 
^mn = ^nm 


where the a*s and d*s are the Fourier coefficients of expansion (A3) 
Then, in the stability equations (AI7) and (A27), 


a, ** 0 


and, for i = J, 


= ^i 


The following upper and lower limits to the true symmetrical 
buckling stress coefficient were computed for the square plate, 
with p = q = 3 in. equations (AI7) and (A27) : 


Lower limit k g = lk. 6 k 

Upper limit k = lk. 79 

B 


Thus, the true value of k s must differ by less t han 1 percent 

from the mean of the upper and lower limits . The numerical results for 
the square plate given by Smith (reference 6) and Iguchi (reference 7 ) 
are lk.72 and lk.58, respectively. 
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APPENDIX B 

FOURIER SERIES IN STABILITY ANALYSIS 


The primary condition in the selection of a Fourier expansion to 
be used as a deflection function is that it be a complete set over the 
region in question. The choice of the best series for a given buckling 
mode, however, depends on certain other considerations which will bo 
outlined herein, with emphasis on application of series in the Lagrangian 
multiplier method. 

Series in one variable It can be shown that any arbitrary 
function f(x) in the Interval (0,a) may be expressed as the sum of 
two other functions, one symmetrical and the other anti synanstrl cal about 

Q 

the midpoint x ■ g- . 

If 

f(x) = S(x) + A(x) 

where 

S(x) = |f(x) + |f(a - x) 

and 

A(x) = |f(x) - ^f(a - x) 


it Is true that S(x) = S(a — x) Is symmetrical and A(x) = — A( a - x) is 
antisymms trl cal . 

When the function f(x) is to be represented by Fourier series, 
both S(x) and A(x) can be given by a series of either sines or 
cosines, as is shown In figure 1. Thus 

S(x) = ^Ta^Sjx) 
m 


Mx) = bjjjAjgC x) 
m 


and 
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Since f(x) = S(x) + A(x), and both S(x) and A(x) acre capable 
of representation by either sines or cosines, there are thus (2 x 2) 
or k possible general series for f(x) that can be constructed from 
the functions of figure 1. 

SerieB in two variables .— Let it be required to represent f(x,y) 
in the region (0,a),(0,b) by a double Fourier series. This double 
series representation can be derived by writing a single Fourier series 
in x and letting the Fourier coefficients depend on y; thus. 


f (*,y) = + t m (y)A m (x) 


m 


m 


Each of the functions and ^ m can in turn be given by Fourier 

series in y, in the region (0,b) as follows 



Then f(x,y) becomes 


fU,y) ■)" a m A m (x)A n (y) Yi W (x)S n (y) 


m n 


m n 


+Y1/L d mn S m (l)S n (y) 


m n 


m n 


It can be proved that 


f (x*y) = S(x,y) + A(x,y) 



to 
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where 


£!(l,7) * a mnV x)A n (y) + 

m n 


51 55 1 »n S » U)S n (y) 

di n 


is symmetrical about the midpoint of the region and 


A( *' y) “5Z2TV A “ (x)S n (y) C mn S m (l)A n ( y ) 


m n 


a n 


is anti symme tri cal about the same point. 

Since the terms S m , S^ A m , in these expressions can be 

given by either sines or cosines from the grohp depicted in figure 1, 
there are thus (2 x 2)^ or 16 combinations of Fourier Beries by which 
a symmetrical or anti symmetrical buckling pattern may be represented. 

For a general pattern, neither symnetrical nor antisymraetrical, a 
choice of 16 x 16 or 256 possible series is available. 

Choice of series .- The simplest application of the Rayleigh-Ri tz 
energy method occurs when the series chosen for the buckled surface 
not only satisfies the boundary conditions term by term, but also leads 
to integrals of products of orthogonal functions in the evaluation of 
the energy expressions. This occurs, for example, when a double sine 
series is used in the analysis of the compressive buckling of a simply 
supported plate (reference 8). 

In more difficult buckling problems, such as c lamped— p la ta problems, 
the simplicity of this calculation can be approached by choosing series 
which do not satisfy the geometric boundary conditions term by term but 
which do have the desirable orthogonal properties ; the Lagrangian 
multiplier method is then used to make the series as a whole satisfy the 
boundary conditions. Thus, in reference 1, a double cosine series was 
employed in the Lagrangian multiplier method of finding .the compressive 
buckling load of a clamped plate. One important consideration to be 
kept in mind in choosing a particular series for use in this n»thod is 
that, for reasons of rapid convergence, the use of cosines rather than 
sines for clamped-edge deflection surfaces is preferable. 

For the case of symmetrical shear buckling of a rectangular clamped 
plate, the following series (equation (A 3 )) was chosen (origin at plate 
midpoint) : 
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W = 





This expansion, having the general form 


^ ^ a mn A m( x ) A n(y) + 
m n 


51 <W B m( x ) s a(y) 

m n 


is a combination of symmetrical cosines and antisymmetrical sines. As 
is desired, the energy expressions involve the products of orthogonal 
functions; the boundary conditions of zero slope and deflection are 
applied to the sines and cosines, respectively, by means of Lagrangian 
multipliers. It was necessary to include sines &b well as cosines in 
the series in order to achieve the desired orthogonal properties; 
however, the portion of the deflection function symbolized 


by 




a S (x)S (y) 
ran m n w ' 


was intuitively believed to be the more 


m n 

important and therefore, for reasons of rapid convergence, was chosen 
to be in terms of cosines rather than sines . Similar considerations 
dictated the choice of equation (A 28 ) as the deflection function for 
antisymmetrical buckling. 
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TABLE 2.- CHOSEN AKRATS 0? SUE AND COS HE TERMS TOR UPPER— LIMIT SOLUTION 
OT SBMETRICAL BUGKLUU, WITH p - 3 AKD 1-2 
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(a) Complete symmetrical sets. 


(0) Complete anti symmetrical sets. 
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